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Abstract 

We study the decay of a very massive closed superstring (i.e. 
a' M"^ ^ 1) in the unique state of maximum angular momentum. This 
is done in flat ten-dimensional spacetime and in the regime of weak 
string coupling, where the dominant decay channel is into two states of 
masses Mi , M2 ■ We find that the lifetime surprisingly grows with the 
first power of the mass M: T = ca' M. We also compute the decay 
rate for each values of Mi, M2. We find that, for large M, the dynam- 
ics selects only special channels of decay: modulo processes which are 
exponentially suppressed, for every decay into a state of given mass 
Ml, the mass M2 of the other state is uniquely determined. 



1 Introduction 



In type II superstring theory in flat, ten-dimensional non-compact spacetime, 
all massive strings are generally expected to be unstable quantum mechani- 
cally by a decay into lighter particles. Massive strings are the key ingredient 
of string theory and crucial for the consistency of the theory. Despite many 
years of study of string theory, very little is known about the way a massive 
string decays. 

Earlier calculations of decay properties are in 0, |^, ^, |, Q , and more 
recent studies can be found in 0, p, ^ . A calculation for states with maximum 
angular momentum in the open string theory was given in An inclusive 
decay rate was computed in 0. A more recent calculation for open and 
closed superstring theory for generic states (which have angular momentum 
much less than the maximum value) is done in 0. 

For small string coupling, the dominant elementary process is the decay 
into two particles. If these particles are massive, then each of them will 
subsequently decay into two ligther particles, and the process ends when 
only massless particles remain. 

Here we shall explicitly compute the rate for the decay of the massive 
string states of maximal angular momentum into two particles. This will 
be done for the closed (type IIA or IIB) superstring theory in fiat ten- 
dimensional spacetime. We will consider all cases: when these particles are 
both massive, when one of them is massless, and when both of them are 
massless. 

A general formula for the decay can be obtained from a one-loop calcu- 
lation of the mass-shift. The inverse lifetime of a massive state of mass 
M is then given by 

ImAM^ (1.1) 

where AM^ represent the radiative correction to the mass, to be expressed 
as a loop expansion. At one loop, it receives a contribution from the two- 
particle intermediate states, and thus T is the lifetime for decaying into two 
particles. 

One can obtain the one-loop expression for AM^ starting from the zero 
and one loop expressions for the four graviton amplitudes. This was derived 
in [§] and we will briefly review this computation in the next Section 2 . We 
get an integral expression in terms of theta functions. [In ^ a computation 
of the lifetime was also attempted; however the algorithm employed in this 
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first investigation was not sufficiently accurate and we find now a different 
and more precise result]. 

One finds that AM^ is formally expressed as a divergent integral of a 
positive quantity: its imaginary part can be computed by analytic continu- 
ation in M^, starting from < where the integral is convergent. By a 
systematic expansion of the integrand the calculation then reduces to a sum 
of integrals of the form 



where each term has a multiplicity depending on a and u. The imaginary 
part is computed by a standard formula: 



One finds that this is precisely the expression for the one loop mass-shift in 0^ 
field theory, and that uj is determined by the masses Mi and M2 of the decay 
products, whereas a is related to their orbital angular momentum, which is 
necessary to compensate the mismatch of J,Ji and J2. This is reviewed in 
Section 3. 

The difficult task is to compute the multiplicity of I{a,uj), due to the 
large multiplicity of the decay products. We have found a very efficient al- 
gorithm for doing that, generalizing the well known saddle-point technique 
for computing the multiplicity (i.e. entropy) of the string states. The gener- 
alities are presented in Section 4. The computation of the imaginary part is 
presented in Section 5. 

In Section 6 we make a numerical analysis of our formulae. This analysis is 
very fast and accurate. We have made in particular a nontrivial check: AM^ 
can be expressed in two formally equivalent ways in terms of two different 
theta functions, which lead to a very different rearrangement of terms in 
the expansion and to different saddle point analyses. We find nevertheless 
identical results, summarized in a three dimensional plot, see figures 1 and 
4. 

We find that Im(AM^) can be written as a sum over the contributions 

of different channels, according to the masses Mi, M2 of the particles of the 
decay product. It turns out that most of these processes are exponentially 
suppressed, except for a line in the space Mi,M2, where the decay rate has 






(1.3) 
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a power-like dependence on the coupling. This is a surprise, since we do not 
see any special reason why, given two masses in a three-state vertex, there 
should be a selection rule for the third. It is likely that that feature results 
from the interplay among the decay products entropy and centrifugal barrier 
and phase space effects Q . 

The other surprise is that we find that the lifetime of the excited state of 
maximal angular momentum for decaying into two particles grows with the 
first power of its mass M. This is computed in Section 7. 

A final remark: at tree level the whole closed string decay can be described 
as a sequence of successive two-body decays. If the decay products are narrow 
resonances - particles with a lifetime much larger than the inverse of their 
mass - then the lifetime of the first decay is the lifetime of the massive state 
tout court, like for instance the lifetime of a i^-meson is set by its decay into 
TT-mesons, despite the vr-meson being also unstable. 



2 The Mass Shift 

Here we review the derivation of the one-loop expression for AM^ given in 
PI by factorizing the one-loop four graviton amplitude obtained by ||10| in 
closed superstring theory (type IIA and type IIB give the same result): 

A, = R' I / n ^ ^-2EL,^.^.x(..) (2.1) 



(Imr)^ J ~}[ Imr 



where 



X{z) = log 



(Im^; 



\2 



2tt^ (2.2) 

Imr 



l^;(0|r) 

and is a kinematical factor containing the graviton polarizations. We also 
set a'M^ = 4N, and choose units where a' = 4, so that in what follows we 
can set = N. 

The amplitude Ai has a double pole for S ^ N {S = 2ki.k2 ) due to the 
propagator of a massive string state with = A^, produced by the collision 
of the two incoming gravitons, and another similar propagator coupled to 
the two outgoing gravitons. The residue of the double pole is proportional to 
AM^, the (1-loop) mass shift of the massive state, ^ + AM^. One 

^We thank Daniele Amati for a discussion on that point. 
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To get AM^, one divides the double pole residue of Ai by the single pole 
residue of Aq, the tree level four graviton amplitude |]TU|: 



An 
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The poles of Ai occur as singularities of the integrand in ( p.lD for zu, z-^i^ 0. 
The integrand behaves as ~ I^i2p'^k34p'^ ■ F (2:12, -234, 2:), with z = + 
^2)- One has to look at the terms in the expansion of F that behave as 
ki2p^~^k34p^~^- Further, in order to select the state of maximal angular 
momentum J = 2N + 2, one looks for the maximal power of cos 6' - the 
angle between the space-momenta of the ingoing and outgoing gravitons in 
the cm. frame - both in the residues of Ai and Aq. In this way one obtains 
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(2.3) 



where s = 7rlm(r) and c is a numerical constant, independent of A^. In 
PI the normalization was checked by evaluating the contribution to ImAM^ 
from the decay rate of the excited state into two massless states, and finding 
agreement with the explicit computation of the decay into two gravitons. 



3 Field Theory 

Let us first consider the case of (j)^ field theory. Consider the one-loop cor- 
rection to the propagator of a particle of mass M, due to the contribu- 
tion of particles of masses Mi, M2 running in the loop. With a convenient 
parametrization, the Feynman diagram has the following form 

AM^ = P ds / dr] e^^^ , (3.1) 

Mi,M2,«o 

with 

W{ri) = M^ri{l -rj)- rjMl - (1 - ri)Ml . 
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Here P{lo) = D/2 — 1 + Iq{D = spacetime dimension), and the polynomial 
P = P{M, Ml, M2, /q) takes into account in particular the multiplicity of the 
decay products, Iq being related to their orbital angular momentum needed 
for matching the angular momentum of the decaying state. 

The IR region is t = cxd. This integral is convergent below the threshold 
for particle production. Above the threshold, the integral is defined as usual 
by analytic continuation, which gives rise to an imaginary part. 

The threshold appears when W{r]) changes sign and becomes positive. 
The maximum of Wit]) is at rj = r]o, with 

_ M2 - Mf + M| 

where 

^(^0) = ^^^'- + ^2)')(M2 - (Ml - M2)') =p\ 

Hence W{r]o) > for > (Mi + Ma)^. 
For future use we define 

+ + )^ .u = ^. (3.2) 

Since AW = ]VPuj — ANPij] — rjo)'^, for large M^, we can evaluate the 
integral over t] by performing a Gaussian integration around the maximum 
at 7] = rjQ. Ignoring constant factors, we get 

'rf^e^^^-^^e^'^^' . (3.3) 

We then evaluate the imaginary part of AM^ by analytic continuation from 
M'^u) < to M'^uj > by means of the formula seen in eqs. ( |1.2[ ) and ( [1.31 ). 
The final result is symmetric in Mi M2. 

In the case of string theory, we will obtain AM^ expressed in two different 
(but equivalent) ways, as discussed in section 4. For comparison with the 
string theory expression studied in detail in appendix A, we set = and 
make the change of variable rj = -y to get 



- iEjpll +47rj/(Af- Af 2 +M| ) 



AM2 = V P r ds s-^^'^^-^e-^^''^^ r dy e 

(3.4) 



Mi,M2,Zo 
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In order to compare with the string theory expression studied in section 5, 
we make another change of variable : 77 = + |, getting the field theory 
expression 

MuM2,lo ° 

" dy e-^-^-^(*^^-*^l) . (3.5) 



2ir 



4 General Method 

Consider the formula for the one- loop string diagram eq. ( |2.3| ), which is 
expressed in terms of 9i{z\t). 

We note that the integrand can be expanded in a sum of terms of the 
form 

where 

can be further expanded in powers of q' = e*^'^'^ and in a Laurent series in 
p = e'™^ , which is symmetric under p —>■ p~^. After the integration over 
Re(r) and Re{z) (ensuring Lq = Lq on the states), we get a sum of terms 
like 

with y = Imz, s = 7rlm(r). 

By comparing with eq. ( p.4|) we see that 2N — (mi + 1712) — 2 corre- 
sponds to k, k to M|, and j to M' - + M| . We can thus determine 
P(M,Mi,M2,/o). 

We will actually compute ImAM^ at fixed Mi^2 summing over Iq, that 
is summing over every possible angular momentum and multiplicity of the 
decay products. 

It is also useful to make a shift z — ^ z + t/2 and to re-express the string 
one-loop diagram in the form 



AM' = c / — ^ / e 



s2 



-mr/4^^(0|r) 
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Now we can expand the integrand of eq. 
form 



(4.4) 

in a sum of terms of the 



T'(mi, 1712) 
where now: 



-2N -2N+(m,i+m2)-3 Ns 



2 .,2 



r 



m2 ' 



Qrr. 
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(4.5) 



(4.6) 



'g-V40i(O|r) 

Note that the new Q'^ has integer power expansions in q and p. By looking 
at the terms \q 
a 
j 



2k 



-2ks 



and IpI"^^ = e '^^^^ (after integrating over Re(r) 
and Re(2;)) and comparing with eq. ( ^.5[ ) we now identify k = + M| and 
Ml - Ml 

The two forms of AM^, eq. (p.3|) and eq. ( [4.4|) are equivalent. However, 
the final expressions that we will obtain using as starting points eq. ( p.3| ) and 
eq. (|4.4| ) and computing integrals by saddle-point approximation will involve 
very different expansions. Therefore it will be a nontrivial check to verify 
that indeed one gets the same result. 



5 Calculation of the imaginary part of 

Here we evaluate ImAM^ using the string loop expression (|4.4|) , written in 
terms of 6*4. 

It should be remembered that the result for AM^ does not change if we 
replace 9^ by 9i, 62 or 6^3, since they differ by a shift of z and a factor that 

compensate the change from e ^ . In appendix A we repeat the analysis 
with AM^ expressed in terms of 6'i, eq. (p.3|) . 

We expand the binomial {^dl log6'4(z|r) + 1)' and using the formula 

^1 1 _ (2iV-mi -ms -2)! 

,=^a.t^„^,) {I - m^)\{l - m2)\{N - I - ir ~ {N-m,-ir{N-m2~ir 

we get (with r = ri + z s/tt and z = x + iy) 



AM^ = c j dsdri j dydx 



V- // X {2N - mi - m2 - 2)\N\'^ , , 

T mi,m2 — — — — T\Ti 
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with T' {nil, 7712) expressed in terms of Qmi2 eq.( |4.5| ). Now we expand 



(5.2) 



with 



Ckj{N, m) 



1 f dq f dp 




(5.3) 



(27r)2 J q J p 



m 1 



and similarly for the complex conjugate. From the explicit expressions for 
the 9i function given in sect. 6, one can see that the sum over k contains only 
positive integer values of fc, whereas the sum over j contains both positive 
and negative powers of j, with the property Ckj = Ck(-j)- 

Let us now consider the integrals over ti and x. Since the imaginary 
part of AM^ comes from the divergence of the integral at s — 00, we can 
replace the integral over the fundamental domain by an integral over the full 
strip. In addition, we note that 64 — ^3 by a shift ti ^ ri + 1. Given that 
the original integral gives the same result for 6^ and 64, we can extend the 
integration region in Ti to the interval (—1, 1). Then 





E e~^^'-^-^y Ckj{N, mi)ckj{N, ms) 



The integration over y is performed by saddle point as in eq. ([3.31). Then we 
consider the integral over s, and use the general formula ( p..3|) for computing 
the imaginary part. We obtain 





V ^' j,fe m,i,m2=U mi,m2=U \ 

Ckj{N,mi) Ckj{N,m2) {Nuj)-^--^ 
mi\{N - mi - l)!2m2!(A^ - ma - 1)!^ 



Here 



uj{p, a) = 1 - 2a + , 



and 



k 




(5.5) 



" N ' 
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We remind that by comparing with field theory we see that the integers k,j 
are related to the masses Mi, M2 of the decay product (cf. eq. 



k = Ml + Ml , i = Ml - Ml . (5.6) 

Here we consider large values of with fixed a and p. Other cases will 
be discussed in appendix B and C. 

It will be clear from the calculation below that the main contribution in 
the sum over mi, m2 comes from the region where 2N — mi — 1712 is large. 
Therefore we can approximate 

(2iV-mi-m2-2)! ^ _ nn _ m2 ^2 

r(2Ar-mi-m2 + |) ^ ^ ' ^ ^ ^ 

Moreover, in the large N limit, the sum over mi 2 is dominated by a sharp 
maximum; away from the maximum the terms are exponentially suppressed 
like e~'^^. We will see that on the maximum r(p, a) = (2 — (mi +m2)/N)~'^^'^ 
is a finite function of j/N,k/N. 
Therefore, we can write 

Im(AM2) ~4-2^Ar2^-i/2^^2iv+| ^^^^^^ \Lij,k)\^ , (5.8) 

where 

(5.9) 

and in turn 

«(9,P,J,*:)- 2. „,(^_„_l),.2 ) ■ (5.10) 



m=0 



6 Numerical evaluation of Im(AM^) 



We now determine the functions appearing in ( ^.8| ). We will make use of the 
formulas in terms of g = e^'^'^^~^ and p = e*^'^^: 



= exp [^/(g,p)] . (6.1) 
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By expanding the logarithms in the above definition of f{q,p), and inter- 
changing the two infinite sums, we obtain 



n=l 

Further: 



fil^P) = -2 E n ^n, ' ■ (6-2) 



(6.3) 

Define 

Now we use the formula ( ^.101) for H = H{q,p,j, k) : 

with n = A^ — 1. For large iV, the integral over t can be computed by a saddle 
point evaluation. We get 

H{y, N) = A^-^+V2 [A^/i(t;)] , (6.5) 

where 

v^m^+i 1 



h{v)=\ogv + \og ^^-—^ ^ + —{VT+4v-l) . (6.6) 

We have checked that this formula provides a very accurate representation 
for the sum in ( |5.10|) already for N larger than 10. 
Thus we finally get 

L{j, k) = ^ ^ ^ ^q-t^p-ie^[fM+hiv(i,p))] (6 7) 

The remaining integrals over q and p can also be computed by a saddle- 
point evaluation. Since the functions appearing in the integrand are com- 
plicated, it is more convenient to perform this calculation by a numerical 
evaluation. 

The saddle-point evaluation of the integrals over q, p is done numerically 
by first finding the extremum of the exponent in eq. ( |6.7| ). This determines 

go = go(p,o-) , Po = Po(p,o-) • 
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We find numerically that the saddle point is obtained for go s'^iid po real and 
positive. By performing the Gaussian integration around the saddle point 
we get 

\L{j, ~ Ar-2^+l"2g27V5,(p,.) ^ (Q_g) 

where 

SLip,cr) = -crlogqo - plogpo + Re[f{qo,po) + h{v{qo,po))] . (6.9) 

It is seen that the saddle is a minimum of Sl for qo,Po real, and that on it / 
and h are real. 

Finally, from eq. (|5.8|) we obtain 



Im(AM2) ~ ^ ^3/2^^^ ^) ^^^^ ^) g2iV5o (g^io) 

where 

S'o(p, (x) = Slip, a) + logt^(p, a) - log4 . (6.11) 



Also, we mention that the result for (|6.6| ) can also be obtained by eval- 
uating the sum over m in eq.( |5.10| ) by looking at the maximum, found for 
mo = A^(l — 2^(Vl + 4f — 1)), and expanding around it (this is done in the 
Appendix C). 

In particular, on the maximum (2 — ^ — ^) = ^{V^ + 4f — 1), since 
the same maximum holds for mi, 2- Therefore we find 

s -7/2 



rip,a)=^-^-—^i^l + Aviqo,Po)-l)j • (6.12) 

Our final formula for the rate of the decay channel to particles of masses 
Ml, M2 is thus given by 



r(Mi,M2) = ImAM^ 
2v 



= ^3/2 1 g2iV5o ^ (5 13) 

Mi,M2 N 



where 5*0 is a function of the ratios with M = \/N. 

Note that r(Mi,M2) represents the contributions of all decay channels 
involving particles with the same masses Mi, M2 (since the multiplicity grows 
exponentially with the mass, for large Mi, M2, there is an exponentially large 
number of particles contributing to r(Mi,M2)). 
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Figure 1 is a numerical plot of 5*0 in function of Mi and M2. We see that 
5*0 is negative definite, except on some curve (see fig. 3) where it identically 
vanishes. Thus the first observation is that the rate for the decay channel to 
particles of generic masses Mi, M2 is exponentially suppressed at large N. 

In the Appendix A, figure 4 shows the numerical plot of So, obtained 
by starting with the expression ( |2.3| ) in terms of 61. The two figures are 
identical. 




Figure 1: The exponent Sq of the deca y rat e r(Mi, M2) in terms of the masses M^/\/N, 
M2/VN, computed using the formula ( |6.11 ). The kinematically allowed region is inside 



the triangle defined by Mi, M2 > 0, ^ +^ < 1- The maximum of 6*0 is iSq = and it 
is located on a curve shown in figure 3. 

Figure 2 is a plot of 5*0 (Mi) for given M2, i.e. slices of figure 1 at constant 
M2. One can see that the maximum exactly passes by Sq = 0. This happens 
for any M2. We have checked that the factor c<j^/^(p, a) r{p,a) in eqs.( |6.10D 
and ( |6.13 ) is finite inside the allowed triangle, except on the boundary Mi + 



M2 = M where F is anyhow suppressed. 

Modulo the exponentially suppressed processes, a massive particle will 
decay through the special channel where 5*0 vanishes. This defines a curve 
M2 = F{Mi) in the space Mi, M2, which is shown in figure 3. Such dominant 
channels exhibit a power-like behavior 

F(Mi,M2) ~ A^-^ . (6.14) 

It is remarkable that for large A^ the dynamics "excludes" decays into kine- 
matically allowed channels. In other words, we find that if the massive 
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particle decays into a particle of mass Mi, the mass of the other particle M2 

is uniquely determined, modulo exponentially suppressed processes. 




-0.2 



-0 . 4 



-0 . 6 



-0.8 - 




Figure 2: Sections of figure 1 at constant A/2- The three curves displayed corespond to 
^ = 0.1,0.3,0.7 (with maxima located at ^ ~ 0.74, 0.28, 0.12 respectively). One can 
see that the maximum always passes by S'o = for every value of M2. 

The curve M2 = F{Mi) is well approximated by the curve 



.Ml 
•If 



M2, 
M 



0.73 



(6.15) 



also shown in fig. 3. Although this not the true analytical formula connecting 
Ml, M2 (which is extremely complicated), eq. ( |6.15| ) is useful as a book- 
keeping of the approximate relation between Mi and M2. 

In the Appendices B and C we have considered the cases when one of 
the masses Mi or M2, or both, are small with respect to A^. We find that 
the only case which is not exponentially suppressed is when M2 = and 
Mf = N — j with j finite (or viceversa) . In this case 

r ~ A^-^/2 _ 



Thus the dominant channel is for both MI2 of order A^ along the curve 



M2 = F{Mi), where the decay width is given by eq. (|6.14| ). 
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0.2 0.4 0.6 0.8 1 



Figure 3: The curve defined by So{Mi,M2) = (solid line), describing the dominant 
channels. For these special values of Mi, M2 the rate is not exponentially suppressed. The 
dashed hne is the curve (Mi/M)"* + {M2/MY = 1, a = 0.73. 

7 Lifetime of a state with maximum angular 
momentum 

Writing the formula (|6.10|) for Im(AM^) in terms of ai^2 = M12/N, we obtain 

Im(AM2) = j^da.da^^^ e^^^" , (7.1) 

where the domain of integration is the region ai, (72 > and ^/Ol + ^Ja^ < 1. 

We have seen that 5*0 is in general negative except on the curve where 5*0 
vanishes. Thus only a small neighborhood of this curve contributes to the 
integral ( [7. 1| ) . Let / be a parameter along the line, < / < 1, and let n a 
parameter for the orthogonal direction, where n = means a point on the 
curve. It is convenient to use n, I as integration variables. The integration 
over / is trivial, since Sq takes the same value (i.e. equal to zero) for all /. In 
the vicinity of the line, we can expand 5*0 in powers of n, and get a Gaussian 
integral over n of the form 

Im(AM2) =N' jdn^^ e"^^"' , (7.2) 
where c is a number of order 1. The Gaussian integral gives an extra factor 
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so we get (see eq. ( |1.1| ) 



= — ^ ^ ^ = const. ^ 



(7.3) 

2\/iV ViV 



or 



r = const a'M . (7.4) 

where we have restored a'. Thus the hfetime of a state with maximum 
angular momentum in closed superstring theory is proportional to the mass. 
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Appendix A Alternative calculation of Im(AM' 

Here we provide an alternative calculation of Im(AM^) by using as starting 
point the formula (|2.3|) in terms of 6i. 

The computation is quite similar to the one done in Sects. 5 and 6, this 
time using the expansion of eqs.( |4.1| ) and ( |4.2| ). We also make use of the 
formulas 



Sm TTZ 



^i(0|r) i\ (1-g 



-,2n\2 



exp[^/(g,p)] (A.l) 



= ~giq,p) . (A.2) 

By expanding the logarithms in the above definition of f{q,p), and inter- 
changing the two infinite sums, we now obtain 

/(g, P) = log[T(p + P-' - 2)] - 2 E ^ n 2n, • (A-3) 
4 ^1 n(l - g2n) 
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Note that / can be defined modulo in, that is modulo the sign inside the 
logarithm, since in eq.( [4.4|) only exp (2Re/) appears. 



Both / and g are even functions of q. We get the same formulas ( |5.9| ) , 
( |5.1CI| ), ( |6.5|) and with f ^ f and g —>■ g and thus v ^ v. Moreover, as 

explained in Sect. 4, the relation with the masses Mi^2 of the exponents j, k 
of the expansions g^^, p' is now different. In terms of the exponents j, k of 
Sect. 5 we have 

~k=ik-j)/2, j = M^-j . 
In conclusion we obtain 

Im(AM2) ~ ^a;='/2(p, a) f (p, a) e^^^^^^''^) , (A.4) 

where 

So{p,a) = -(cr-p)loggo- (1 -p)logpo 

+Re[/(go,Po) + h{viqo,Po))] + log t^(p,cx) , (A.5) 

_ + M| _ Mf - M| 

and the saddle point values qo,Po correspond to the stationary point of Sq. 



The functions that enter into the final expression (|A.5|) are individually 
very different from the ones appearing in the formula for 5*0 in Sect. 6. The 
result is however the same, to a surprising degree of accuracy, see fig. 4. 



Appendix B Decay into two massless particles 

The decay rate of the massive particle into two massless particles (e.g. gravi- 
tons) can also be obtained from the general formula ([4.4|) . We need to con- 
sider the special channel with Mi = M2 = 0, i.e. the powers q'^p' in the 
expansion with k = j = 0. This is obtained by formally setting g{q,p) — > 
(since it starts with q^) and ^-iftg^jj^^^^.^ ~^ 1- For this k = j = term, we 



thus get 



massless 



dsdn f <fz ^ 1 N\2N-2)\ 



J s 'As' {N-iy/ 
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Figure 4: The exponent Sq of the de cay r ate r(Mi, M2) in terms of the masses Mi/M, 
M2/M, computed using the formula ( ]A.5| ). The plot should be compared to figure 1, 
computed by following a different way. 



Computing the integrals over x,y and ti, we are left with an integral over s 
whose imaginary part gives 

Im / ds s-t-2^e^^ ^ ^e2^2-2^ 



Hence 



ImAM^ 



~ ^/N e-2^('°*5^-^) . (B.2) 

massless 



This is exponentially decreasing for large A^. 

Note that this decay corresponds to the corner Mi = M2 = in figure 1, 
which is also exponentially decreasing. Indeed, the numerical value at Mi = 
M2 = of the plot is = -0.7726, which agrees with -2(log4 - 1). 

In § this result has been compared with the explicit direct computation 
of the decay into two gravitons, finding complete agreement. 



Appendix C Decay into a massless and a massive 

particle, and remaining cases. 

Another interesting special case is when the massive particle decays into a 
massless particle (e.g. graviton) and a massive particle. So, consider the 
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decay M ^ Mi + M2 with M2 = and = N - j (we remind that 
= A^). We consider both the case when j is a finite integer and the case 
j = c N where c can be a constant less than one or c ~ iV'^~^ with A < 1. 
The case M2 = is most suitably treated by the formula ( |2.3| ) expressed in 
terms of 6i{z\t), setting g = to isolate the term corresponding to M2 = 
in the expression ( [4 .31 ). Thus (see (|A.1|) , ( |A.2| ) ) 



sin(7r2;) , —d^\og6i{z\T) 



9[{0\t) ' ' ' 7r2^ ^ ' ' sin2(7rz) ' 

Therefore now 

4 

uip) = g{q = 0,p) = ^^i/2_p-i/2)2 ' 



e 



/(9=0,p) _ l^M^X 



e'Mr) 



-u{p)-' . (C.2) 



Now u = thus ujN = j'^/N and ^ = ^^M^. By following similar steps 
as those which led to eq.(p).4[), now we obtain 



IniAM^ ~ M)L^/^^p-^(«(p))-^/^p-^(«(p))-^ 
\/N J V J V 



N J p J p 

^ X! rToM ~ ^ , \\ -^mz • (C3) 

mi ,r7i2 = 



where 



= Z{u{p),m,j,N) 



r(2A^ - mi - m2 + f ) 

m Mp)n^ 



We use a saddle point technique to sum over mi ^: the dominant contribution 
comes from the maximum of the exponential dependence in those variables. 
Using the Stirling formula, we find 

Z^r^{N-m-l)M e^i°siVe^ , (C.4) 
V m 

I = N -m-2Nlog{N -m-1) 

+ m f-logm + 2 log (iV-m- 1) + log^^^^^'j . (C.5) 
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Imposing that the derivative in m of the exponent vanishes, we get the equa- 
tion for mo, the maximum locus. 



(AT - mo - 1) 



,2 ^Oj 



u{p)N ' 



which is solved to give 



Let us now discuss different cases. 
For f/N ^0 , eq. ( [(TBI ) gives 

mo = - ^ - 1 . (C.7) 

On the maximum, we get 

Z{u{p),mo,j,N) ~ -^{^f «^ exp[-^] , 



{2N -nil - ma -2)! .^7/2/ 1 ^ 1 y7/2 
r(2iV-mi-m2 + |) 



and 



The second derivative of the exponent / is — 1/mo — 2/(A^ — mo — 1), which in 
this limit is ^yfuj j. We are left with a Gaussian integral in 5mi, bm2 with 
a spread of order 1/7, which is small compared to the range of m. Hence 



ImAM^ ~ v-^v (-^27v -7/2 



J ' 9 — p exp —= 

where we have neglected finite powers of u and kept into account the range 
m m. We recall that 1/01 = _ p-i/^y xhe remaining integral over 

p is obtained again by saddle point technique: defining x = p^^^, we require 
the derivative in x of the exponent to vanish, that is —2x^^ + |(1 + x~'^) = 0. 
The solution is a; = Xq = 2 + y^, discarding the other solution x = I/xq 
which would make N — m ^ j/u{p) < 0. We get 

f exp [^{p^^^ - p'^^^)] ~ exp [ - coj] , 
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where 



Co = -21ogxo + ^2^ = 0.9 . 

ZXn 



Thus we finally obtain 

•3/2 

ImAM^ ^ ^ exp [ - 2coi] • (C.9) 

For large j, this rate is exponentially suppressed. 

The calculation applies as well to the case of finite j. In this case, the 
exponential factor depending on j is a finite number of order 0(1), so one 
obtains 

r ~ _ 

Now consider the case when N > j > N^'^. Then eq. (|CJ) gives N — 

■1 

mo — 1 ~ ^ — and therefore we find the same result, up to negligible 
corrections. 

Finally, consider the case j = cN. Then eq. (|C.6| ) reduces to 



1 Nc^ / 4-u 

iV-m=^(^-2)(Jl + ^-l) 

and we have to evaluate 



P 



J — p exp 



-logu-21o£ 



1 + f - ly 

2u 




We note that this is a particular case of the expression of the general case 
(|6.71 ) - we recognize v{q = 0,p) = u{p)/c^ and the expression N{f{q = 
0,p) + h{y{q = 0,p)) in the exponent - except that in the general case 
M| ~ iV we have an extra Gaussian integral coming from the integration over 
q (rather than setting g = to isolate the term corresponding to M| = 0). 
As a result, ( |6.7| ) will give an extra factor |l/v^p. 

Summarizing, for M2 = and Mf = N — j with 1 < j < cN, we get 

■3/2 

ImAM^ ~ ^ exp [ - bj] (C.IO) 

where 6 > is a numerical coefficient of order 1. 

The case when j = N, i.e. c = 1, corresponds to the Mi = M2 = 
case treated in appendix B. It can also be recovered as follows. We need to 
compute the coefficient of the power p'^ in the Laurent expansion of 
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exp [N{ - logM - 2 log ^/"^ ^ + ^{VTTi^ - 1))] 

= (V exp [iV(l + ^ c„p-")] (C.ll) 

n>0 

The coefficient of is (|)^. Now we have to take the modulus square of 
this coefficient, and take into account that the generic case eq.(|C10[) has 



an extra factor originating from the Gaussian integral around the 

saddle point over p. We thus obtain the correct result for Mi = M2 = 0: 

ImAM^ ~ . , (C.12) 

in agreement with appendix B. 

Finally, we have also considered the case when M| = n is finite or small 
with respect to N, which is allowed for j > — n. We find that this 

case is exponentially suppressed. 
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